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ABSTRACT: It is shown that some formulae describing optical rotatory dis- 
persion of a-quartz with the aid of two Drude's terms reduce to the combined 
formula containing one Drude's and one Chandrasekhar's term. Comparison of 
various formulae describing the experimental data of a-quartz leads to the con- 
clusion that the optical activity of this crystal is due to its crystal structure only, 
that means the optical activity is not of molecular origin. Further the rotatory 
strengths are discussed with the regard to coupled oscillator model and to the 
structure of a-quartz. 

1 Introduction 

Burkov et alii jl], |2j have undertaken an exhaustive study of formulae describ- 
ing optical rotatory dispersion (ORD) of a-quartz. This study is based on precise 
measurements of various authors || |4], [5|] in the region from 0.15/im to 3.2/xm. 
Beside the well known formula of Drude 

K {1) 

Burkov et alii have used the following formulae: 

^ 1} A 2 
(A 2 - A 2 ) 

derived by Chandrasekhar || 0, the combined formula 



p(A) = Ett^W ^ 



P(A) = E 



+ 



A 2 - A 2 (A 2 - A 2 ) 2 



(3) 



1 



derived by one of us || and also by Nelson |J and formula 



P(A) = E 



iff (A 2 + A 2 
(A 2 - A 2 ) 2 



(4) 



derived by Agranovich 



and by Tsvirko ||12|| . In these formulae p is the 
rotatory power, Aj is the characteristic wavelength and are the constants 
including molecular and crystalline properties of a-quartz. Formulae (||]) and 
(0) are based on the theory of coupled oscillators and the formula (|J) is based 
on the excitons theory. It should be noted that measurements and theories for 
light propagating along the optic axis have been taken into account only. The 
applicability of the above mentioned formulae have been judged by Burkov et alii 
with the aid of root-mean-square deviation defined by 



9 ^} j yPexy Pteor) 



(5) 



where p exp are experimental data and pteor fitted data following from formulae 
(0) " (i)- Further n is a number of measured data and the analysis of Burkov et 
alii comprises forty-five points from 0.15pm to 3.2pm. These authors have found 
that formulae ([!]) and (Q) with one term are inconvenient because g is very great. 
They have found good agreement in the case of formula (|I|) with two terms of 
opposite signs, in the case of formula (H) with one or two terms and in the case 
of formula (§) with different and K% but with one Aj. In some cases they 
have introduced the constant term corresponding to ORD in infra-red region. 
The results of the most convenient formulae are presented in Tab. [I} 



Table 1: Table of formulae. 



Nr. 



Formula 



7.17A 2 ~~~ 
(A 2 -0.0928 2 ) 2 
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I 

II 
III 

IV 
V 
VI 



P(A) 
P(A) 
P(A) 
P(A) 
P(A) 
P(A) 



-0.05 



7.266A 2 



A 2 -0.09266 2 

184.22 
A 2 -0.09355 2 

-1.05A 2 



+ 

+ A 2 -0.09170 2 



(A 2 -0.09266 2 r 
177.04 



I 8.22A 2 
(A 2 -0.09998 2 ) 7 (A 2 -0.09387 2 ) 2 



219.47 



A 2 -0.09332 2 
-1.335A 2 



+ 



-212.30 



(A 2 -0.1012 2 ) 2 ~"~ (A 2 -0.0944 2 ) 



A 2 -0.09176 2 
8.49A 2 



+ 0.0361 
+ 0.0964 



3.91 
3.90 
3.90 
3.82 
3.89 
3.74 



In addition it should be noted that Katzin and Burer |T3| , [14| have worked 
out a formula for ORD of a-quartz in the range from 0.23pm to 3.5pm. Their 



2 



formula has the form 

P(A) = 



127.02476 
A 2 - 0.0979 2 



(6) 



The experimental data are fitted to a root-mean-square deviation of about 
0.038deg. It is the most accurate formula worked out for this spectral range. 
However, this formula deviates notably from important measurements of Servant 
|| in the far ultraviolet. In the range from 0.1850/im to 0.15235/im the deviations 
are of systematic nature and become very great. E. g. for A = 0.15235/im rotatory 
power is 779.9deg/mm whereas formula (6) gives 791. 7 deg /mm. 

The good agreement of formulae containing the wave dependence A 2 / (A 2 — A 2 ) 2 
is not surprising because this dependence holds for some other crystals which are, 
like a-quartz, optically active in crystalline state only. On the other hand good 
approximation with the aid of formulae containing Drude's terms only such as the 
formulae III and V from Tab. [l] is surprising because we are suspicious, in agree- 
ment with Hennessey and Vedam [Tj|, of the formulae containing Drude's term 
only. These terms are typical for optically active molecules. The molecular origin 
of optical activity of a-quartz is highly improbable because of the high symmetry 
of S1O2 groups. Of course, they are some crystals the ORD of which can be 
fitted by Drude's formula although they consist from symmetrical molecules. In 
this case the optical activity is of a pseudomolecular origin. The molecules of 
these crystals are symmetrical when isolated but they lose their symmetry in the 
crystalline state. They play then the same role as chromophors in optically active 
molecules. But it seems that this situation occurs in the case when crystals are 
composed from large molecules. We shall see that the loss of symmetry of SiOi 
groups in a-quartz is very little. We shall show that the formulae containing only 
Drude's terms reduce to our formula (^), where Chandrasekhar's terms plays a 
predominant role and Drude's terms represents a little contribution. It means 
that formulae III and V differs very little from formula II. 



2 Coupled Oscillators Model 

ORD of a-quartz via coupled oscillators model has been firstly studied by 
Chandrasekhar || [7j and then by us |§. This model has been also generalized 



to the case of absorption region, where circular dichroism occurs |16], [P7|| . In 
this section we shall improve this model in order to obtain equations containing 
explicitly the rotational strengths for normal modes of vibrations which are the 
specific features of this model. We shall restrict our considerations to the crystals 
structure belonging to the space group of symmetry D| or its enantiomer 
because a-quartz is its typical representative. The crystal structure of a-quartz 



19] shows that an axial channel which results in a lower specific gravity than 



would be found from closer packing of the atoms. Fusion, rather than giving a 
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higher density from closer packing, results in an even lower specific gravity. This 
implies a molecular S1O2 structural unit rather than an ionic crystalline constitu- 



tion [18. 19]. Unit cell contains three molecules which are arranged spirally about 
a set of parallel c-axes. This arrangement together with the interaction between 
molecules causes a gyratory effect. In order to describe the optical properties 
of medium under the consideration it is natural to start from the properties of 
an isolated molecule and interactions between molecules are neglected because 
they play the role of little contributions. But it is not possible in the study of 
the optical activity of crystals consisting from inactive molecules. Therefore the 
above mentioned chiral structure together with the interaction between molecules 
must be simultaneously taken into account. In the study of interaction of an iso- 
lated molecule with an electromagnetic radiation each molecule behaves like a 
linear harmonic oscillator by which it may be represented. In the study of the 
optical activity the interaction between molecules can be considered with the aid 
of a compound oscillator consisting from two linear harmonic oscillators coupled 
together. 

According to crystal structure of a-quartz let the first single oscillator, rep- 
resenting an isolated Si0 2 group, be situated at the origin of co-ordinate system 
with its vibration direction along a, /3, 7, where a, /3, 7 are the direction cosines. 
Each oscillator is a charged particle which is assumed to be bound elastically to 
its own equilibrium position and capable of vibrating along a line. In the case of 
a-quartz we assume the vibrations along the line joining the centers of gravity 
of oxygen atoms. The second oscillator is situated at x, y, z with its vibration 
direction turned through an angel 9 about c-axis with respect to the first one. For 
a-quartz 9 = 120deg. Each oscillator, when uncoupled, has the natural frequency 
ujq. As a result of the coupling, u would be split into two normal frequencies uj\ 
and o>2- These normal frequencies of the compound oscillator are 

uj\ = + 27r 2 e, = Uq — 27r 2 e. (7) 

Here e is a coupling constant. The exact nature of the coupling is not essential 
in the following considerations but we consider it as small. We can neglect all 
terms containing e 2 . In our case it means that we assume the interactions of 
valence shall electrons of each isolated Si0 2 group. 

The above mentioned model enables us to evaluate ni and n r as it has been 
performed by Chandrasekhar H [7| • He has assumed that both normal modes have 
the same oscillator strengths. We have removed this simplification || obtaining 
for rotatory power p the following expression 



nNe 2 {a 2 + (3 2 )dsin9 ( f q2 uj 2 f qi u 



,2 



me- \ ui 2 — us 2 io\ — uj 2 



(8) 
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or in wave lengths 



P(A) 



nNe 2 {a 2 + f3 2 )d sin 



A 2 - A 2 



A 2 -A 2 , 



(9) 



From eqs. (0) and (H) one gets 



P(A) 



TiNe 2 (a 2 + (3 2 )dsin6 



e(U + U)XtX 2 , (/, 



'/2 



L 2c 2 (A 2 -A 2 ) 



2^2 



A 2 -A 2 



(10) 



Here N is a number of single oscillators per unit volume, d is a linear distance 
between two single oscillators and f qi , f q2 are oscillator strengths of both normal 
modes of vibrations. 

We see that this model does not contain rotational strengths explicitly. It 
is because this model is a classical one. The rotational strengths follows from 
the quantum mechanical equations of motion of coupled oscillators. It will be 
shown in Appendix. Each normal mode has a proper rotational strengths. The 
strengths are 

he 2 d(a 2 + (3 2 ) f qi sm9 



R 



91 



R 



8mc 

he 2 d{a 2 + [3 2 )f q2 sm6 



1-12 



Substituting from eq. (11) into eq. 



8mc 

3) we have 



P(A) 

or with the aid of eq. (0) 



8vriV 

he 



Rq 2 ^2 

A 2 -A 2 



+ 



Rqi^l 

A 2 -A 2 



P(A) 



AnNe (R q2 -R qi )X 



X 2 



he 3 

8nN (R qi +R q2 )X 



(A 2 - A^ 



he 



A 2 -A 2 



+ 



2 ' 



(12) 



(13) 



It is clear that eqs. fljjD and @ are identical but under the condition that 
the difference between Ai and A2 is very small. This little difference follows from 
the assumption that e is small. Now we see that eq. ( |13"1) contains two terms 
the first of which is of Chandrasekhar's type while the second one is of Drude's 
type. In the case of f qi = f Q2 we have R qi = —R q2 and Drude's term in eq. (pf) 
vanishes. Comparing eqs. ( |T2"D and ( |T3| ) with the formulae in Tab. [I] we see that 
formulae III and V, being of Drude's type, contain Ai and A2 which are very near 
and that the constants are also very near in numerical values but of opposite 
signs. We can also expect that these formulae reduce to our combined formula 
(Pf). Let us suppose that both constants in formulae of the type of III and V are 
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the same in numerical values but with opposite signs, then these formulae reduce 
to Chandrasekhar's one. We can also expect that formulae III and V reduce to 
the formula which should be very near to formula II of Tab. [l|. We shall see in 
the next section that this is true. 



3 Discussion of Formulae 

Let us compare formula III with eq. (|12|). We rewrite eq. ( |12|) as follows 



We must take care of units. In formulae I - IV p is expressed in deg/mm and 
in theoretical considerations we express it in rad/cm. In eq. ([14]) p and A q2 and 
A qi are in the same units. Comparing eqs. (|l^) and (|i~4"l) we see that 



91,2 = -sr 4 - < 15 > 

Substituting from eq. ( |T5D for R qi — R q2 and for R qi + R q2 into eq. ( p~3| ) one 

gets 

(eA 4 /2c 2 ) (A q2 - A qi ) A 2 , (A qi + A q2 ) A 2 

P(A) " (A^T^ + A 2 - Aq ' (16) 

where A = (Ai + A 2 )/2. Using eq. © we get 



i 



eA 



I (A 2 - A 2 ) = A AA. (17) 



2c 2 

Here AA = A 2 - X v It follows from formula III that A q2 = 2.1049 ■ 10 4 and 
A qi = -2.1053 ■ 10 4 in units of p. Further A = 0.09265/im and AA = 1.85nm. 
Substituting these values into eq. (|I6|) we have 

7.214A 2 0.034 . , 

p(\) = ? - — (18) 

PV ; (A 2 - 0.09265 2 ) A 2 - 0.09265 2 V ; 

In a similar way we obtain instead of formula V 

,(A)= 728U2 °' 013 2 + 0.036. (19) 

HK ' (A 2 - 0.09254 2 ) 2 A 2 - 0.09254 2 v ; 

In this case Ao = 0.09254/im and AA = 1.54nm. We are able to evaluate 
e from eq. (|T7D . It is 4.25 • 10 29 sec~ 2 for formula III and 3.54 • 10 29 sec~ 2 for 
formula V. Because Uq is of order 10 32 sec~ 2 we see that e is really small as we 
have assumed. Chandrasekhar H has also evaluated e for a-quartz but he had to 
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use beside his formula I the formula for ordinary refractive dispersion. His value 
is 6.35 • 10 29 secr 2 . 

We see that there is a little difference among formulae I, II, III and V. We 
conclude that formulae III and V are nothing but formula II. 

Now let us discuss the rotational strengths which may be evaluated from eqs. 
( |T2"D and (0). In this case we must multiply both sides of eq. (|14|) by 7r/18. 
Then we obtain p directly in rad/cm, that is 

3.664 ■ 10 3 A1 3.674 ■ 10 3 A 2 
= A 2_ A 2 A 2_ A 2 ( 20 ) 

or 

3.674 -W 3 hc 3.664 • 10 3 ftc 

^ " 8^iV~' R< * ~ 8nN ■ (21) 

We have for a-quartz N = 2.68 • 10 22 and therefore R qi = —1.732 • 10~ 37 and 
R q2 = 1.727 • 1(T 37 in cgs units. 

We are now able to evaluate oscillator strengths f qi and f q2 . We have for 
a-quartz § a 2 = 0.64, (3 2 = 0.0232, d = 1.8 • 10~ 28 cm, sm6 = y/3/2. We get 
from eq. (|TT| ) 

R qi = -1.155 • 10" 38 / gi = -1-732 ■ 10" 37 
R q2 = 1.155 • 10" 38 / g2 = 1-727 • 10~ 37 . (22) 

From thus f qi = 14.99 and f q2 = 14.95 and the average value fo = (f qi + 
f q2 )/2 = 14.97. This actual value of oscillator strengths corresponds to the num- 
ber of valence - shell electrons of the oxygen and silica bring into S1O2 molecule. 

Using the same procedure in the case of formula V we get for R qi = —2.0675 • 
10~ 37 and for R q2 = 2.0677 • 10~ 37 . The splitting of oscillator strengths is very 
little because f qi = 18.852 and f q2 = 18.851. For these reasons Drude's term plays 
a negligible role in formula V. The averaged value of fo = 18.8515 seems to be in 
high regard to the number of valence - shell electrons in S2O2. Chandrasekhar 
combining formulae for ORD and refractive dispersion has obtained 

^~Y^h = 6-5. (23) 

From thus fo = 19.68 which is also higher then the valence-shell electrons in 
SiQ 2 . 



4 Conclusions 

Performing the analysis we may conclude that all formulae from Tab. [l] reduce 
to formula in which the term A 2 /(A 2 — Ag) 2 plays a predominant role. This 
wave dependence arises from the model of coupled oscillators. Although this 
model is not sufficiently general to be rigorously valid for a crystal, it brings out 
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clearly a type of interactions which may be important in determining ORD of 
crystals which are optically active in crystalline state only. It has been shown 
that formulae containing two Drude's terms with opposite signs and involving Ai 
and A2 which are close to each other reduce in fact to combined formula in which 
Chandrasekhar's term plays a predominant role. 

Nelson || has analysed the mechanism of crystalline optical activity and con- 
cluded that the quadratic term in formula II is due to lattice dynamics, while the 
electric-dipole-magnetic-dipole or electric-dipole-electric-quadrupole interference 
leads to the known Drude's term. 

The evaluated rotational strengths of both modes of vibrations are of an order 
10~ 37 . This value is acceptable for crystals in visible and in uv region. On the 
other hand it seems that only the formula III and for these reasons also the 
formula II gives the actual value of /o corresponding to the number of valence - 
shell electrons in Si0 2 - 

The location of Ao is also of great importance and facilitates to understate 
the mechanism of ORD in a-quartz. Platzoder has measured the vacuum uv 
reflectance of a-quartz and has found that the sharp peak at 0.1190/im is due 
to this exciton transition and that the other peaks at 0.105/xm, 0.0840/wn and 
0.0710/im are due to the interband transitions. The characteristic wave lengths in 
formulae I - III and V are located at the bottom of the dip between two peaks at 
0.084 and 0.105/im. The bottom is at 0.0940/im and it is in agreement with ORD 
formulae. On the other hand characteristic wave lengths in two - termed formulae 
IV and VI belong to two different interband transitions. The peak at 0.105/im is, 
in fact, split into three very near peaks. Two bottoms are located between 0.10 
and 0.11/zm. The characteristic wave lengths in first terms in formulae IV and 
VI are then very close to above mentioned values. Both other terms contain Ao 
which is located at « 0.084/im. 

The connection between characteristic wave lengths and interband transitions 
is in agreement with the model theory of coupled oscillators. The normal modes 
of vibrations are nothing but the tight-binding exciton model in quantum me- 
chanical term. 



Natori JH] has pointed out that coupled oscillator model pays no attention 
to the effect of the electron transfer from atom to atom that gives rise to the 
ordinary electronic band. For this reason Natori has discussed ORD in term of 
the ordinary energy band theory. He has obtained for a-quartz following formula 

p(u) = 911.22(^10.834 - Tiuj + VlO.834 + Tiuj - 2^10.834) + 

1 12 

+5708.8 f ; = + ; ) , (24) 

Vl0.834-faj V10.834 + Tiuj VlO.834 7 

where Tiuj is measured in eV. In this formula he has taken into account 11 points 
of measurements from 0.6708 to 0.1525/zm and calculated p with the aid of eq. 
(|24|). But in this region g = 2.16 while Chandrasekhar's formula gives g = 0.992. 
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Moreover it follows that the deviations Ap = p exp — p ca i are of systematic nature 
when using formula (p4|) . In the region from 0.6708 to 0.3403pm Ap in negative 
and the deviations have decreasing tendency from —0.11 to zero at A = 0.3403pm. 
In the region from 0.3403 to 0.1525pm Ap is positive with increasing tendency 
from zero to +3.9 at A = 0.1725pm. The deviation Ap at A = 0.1525pm is again 
negative but very great, that is —5.8. In formula (25) Hujq = 10.834el^. It means 
that Ao = 0.11443pm which located it at the bottom of the dip between two 
peaks at 0.105 and 0.1190pm. Natori has pointed out that this location of Ao is 
in agreement with the used model and that the formula (|24|) is more preferable 
than Chandrasekhar's formula. But as we have shown his formula does not enable 
good approximation of experimental data. 

The analysis cannot stop here but we can conclude that ORD of a-quartz is 
due to its crystal structure because in all appropriate formulae Chandrasekhar's 
terms play a predominant role and just these terms reflects the optical activity 
of crystalline origin. It should be noted that Drude's term in our formula does 
not reflect optical activity of molecular origin as judged by some authors (see e. 
g. [2^] and references in this paper). This term reflects the splitting of rotatory 
strengths of normal modes of vibrations (see Appendix). Meanwhile, we do not 
know the cause of this splitting but we expect that this splitting is caused by 
crystalline field. We hope in the future to apply the model of coupled oscillators 
to some real crystals including the effects of crystalline field. 

Some authors p2|| ascribe our combined formula to Agranovich [ 10|, |i~l~| . But 



Agranovich has pointed out that Drude's term vanishes in his formula when 
taking into account crystal consisting from optically inactive molecules. The 
difference between our formula and formula of Born j33| has been discussed by us 
in P]. In the Appendix we shall show that formula of Agranovich and Tsvirko is 
special case of our combined formula as well as formula derived recently by Kato 
et alii p||. 



Appendix: Rotational Strengths of Normal Modes 
of Vibrations 

A transition is optically active when it has non-vanishing rotational strengths. 
For a transition from some ground state |0) to the excited state \n) the rotational 
strength is a imaginary part of scalar product of the electric dipole moment p e 
and the magnetic dipole moment p m . These moment are (0|p e |n) and (n|p m |0) 
so that 

Ro n = lm(0\p e \n)(n\p m \0). (25) 

In the case of coupled oscillators the excited state |n) is split into two states 
| In) and |2n). Therefore the rotational strength needs be considered between 
ground state and both excited states. The total electric moment of compound 
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oscillator is 



fie = e (ft + f 2 ) = 

= e[mri + jflrx + i(a cos 6* — (3 sin #) r 2 + J(a sin 6 1 + (3 cos #) r 2 ]. (26) 

Here ?and J are the unit vectors along x and y axis and r*i and r 2 are position 
vectors of vibrating electrons in both single oscillators. Because we study the 
rotational strengths for light propagating along c-axis we taken into account the 
vibrations along x and y axes. Using the normal coordinates 



ri = (<?i + <?2) , 



we obtain from eq. (|26|) the normal components of follows 

j2 £qi = ~^{^[ a (1 + cos ^) ~~ P s i n 0] + (1 + cos#) + a sin 0] jgi 



and 



/L = -^={t[a (1 — cos6>) + {3 sin 6 1 ] + (1 — cos^) — asin6>] } 
2 V2 L J 

The magnetic dipole moment of compound oscillator is 



Q2- 



(27) 



(21 



(29) 



2c 



x r 2 



x ri 



(30) 



where <i is linear distance between two centers of gravity of both single oscillators. 
We can write instead of eq. fl30l) 



ed 

fIm = ~Tc 



k x m — 



r 2 



(31) 



Here k is the unit vector along z-axis which is identical with crystal c-axis. 
We see that magnetic dipole moment arises from a particular spatial distribution 
of coupled oscillators. When f\ and r 2 colinear the magnetic dipole moment is 
zeroth. The normal components of 



and 

— * 



ed 



32c 



[j[a (1 - cos 9) + (3 sin 9] - i\J3 (1 - cos 9) - a sin 9] }qx (32) 



ed 



32c 



[j[a (1 + cos 9) - (3 sin 9} - i[(3 (1 + cos 9) + a sin 9} }g 2 . (33) 
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According to eq. (|25| ) we get 



Ri n0 = Im 



and 



R2 



Im 



2e 2 d(a 2 + (3 2 ) sin0(O|gi|ln)(lra|gi|O) 

8c 



2e 2 d(a 2 + (3 2 ) sin 6($\q 2 \2n) (2n|g 2 |0> 
8c 



The matrix elements {f]n\q v \Q) may be written as follows 

(?H<yO) = ^no(^|^|0) 
and using definition of oscillator strength 

2mu vn0 \(r}n\q 11 \0)\ 2 



we obtain for R ln0 and i? 2 „ 



and 



#2 



he 2 d (a 2 + p 2 ) f ql sin 9 



he 2 d(a 2 + (3 2 )f q2 sing 
8mc 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



We see that R\ n0 and R 2n0 substituted into eq. (|T^) gives directly eq. (|9|) and 
from eq. (p~3| ) we get eq. ([U]). 

In linear harmonic oscillator approximation one may evaluate the matrix el- 
ements in eq. (1371), that is 



\{vn\q v \0)f 



2muo. 



(40) 



and from thus f qi = f q2 = 1. Then Drude's term vanishes in eq. (10) and we get 
Chandrasekhar's formula. Of course, in real molecules f qi = f q2 < n, where n 
is a number of oscillating electrons. In Heitler - London approximation one may 
write (ln|gi|0) = (2n\q2\0) and therefore 



fqi _ fi 



92 



A 



(41) 



Here /o is> as before, average value of both oscillator strengths, that is os- 
cillator strength of single oscillator. Further u) n o is characteristic frequency of 
uncoupled oscillator. Using eq. (0) we get from eq. (|4l|) 



fqi + fq2 — 2/o 



fq2 f' 



qi 



2rf 2 f e _ ef X 2 
2c 2 



<^n0 



(42) 



11 



and substituting these expressions into eq. ( |T0| ) one has 



nNe 2 d{a 2 + (3 2 )ef \ 4 sm6 A 2 + A 2 , 



mc 

or 



; -" 4 (A 2 - Aq) 

( 3 ) (A 2 + 
(A 2 - \h 



(A 2 + A, 



2\ 



pw = : 2 ; 2 o; - (44) 



This formula is identical with that obtained by Agranovich [IT], [TTJ and by 
Tsvirko jnj. We see that formula ( |4*4] ) is nothing but the special case of our 
combined formula (f|). 

Kato et alii have derived ORD formula from the exciton dispersion theory. 
Their formula is 

= ww-® 

HK ' (A 2 -A 2 ) 2 1 ; 

Without entering into the details of the microscopic theory it is interesting 
that when = 2K<U in our formula (^) this formula reduces to formula (fl^p. 
On the other hand formula (^) may be formally obtained from formula (|3[) when 
2K( 1 ' = —K^ but it is nothing but Heitler - London approximation discussed 
in eqs. (P) - (M). 
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